for k=1:45
    
    % Todo: 后续找5个样例验证
    P(k)=sum(sum((abs(DIFF)<k)))/(size(DIFF,1)*size(DIFF,2));

%     title('Lateral shooting','Interpreter','Latex')
end
plot(P);
set(gca,'XDir','reverse')
xlabel('K','Interpreter','Tex')  %x轴坐标描述
ylabel('P','Interpreter','Latex') %y轴坐标描述